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Abstract 

We derive the statistical distribution functions for the Hubbard chain with 

infinite Coulomb repulsion among particles and for the statistical spin liquid 

with an arbitrary magnitude of the local interaction in momentum space. 

Haldane's statistical interaction is derived from an exact solution for each of 

the two models. In the case of the Hubbard chain the charge (holon) and 

the spin (spinon) excitations decouple completely and are shown to behave 

statistically as fermions and bosons, respectively. In both cases the statistical 

interaction must contain several components, a rule for the particles with the 

internal symmetry. 
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I. INTRODUCTION 



It is well known that in the space of dimension higher than two the many-particle wave 
function is either symmetric or antisymmetric under a permutation group operation; this 
property leads to the division into the systems of bosons and fermions, respectively. As a 
consequence, the distribution function for the ideal gas is given either by the Bose-Einstein 
or by the Fermi-Dirac functions ||I[]. In low-dimensional systems {d = 1 and 2) the situation 
changes drastically because e.g., a proper symmetry group in two dimensions for the hard- 
core particles is the braid group, the characters of which are complex numbers |0. In such 
instances the distribution function has not been determined as yet. On the other hand, the 
distribution function can be changed by the interaction among particles. Such a situation 
arises for instance at the critical point when the system undergoes a phase transition. Be- 
low the critical temperature (e.g., in the superconducting phase) the distribution function 
changes its form from that in the normal state. So, the statistical properties of the particles 
are influenced by both system dimensionality and by the character of dynamical interaction 
between particles. 

In his remarkable paper 0], Haldane noted that the distribution function can also differ 
from the Bose-Einstein or the Fermi-Dirac form in the normal state. He generalized the Pauli 
exclusion principle by introducing the concept of statistical interaction which determines 
how the number of accessible orbitals changes when particles are added to the system. The 
paper dealt with the many-particle Hilbert space of finite dimension. The limitation turned 
out to be irrelevant. Namely, Murthy and Shankar showed that the statistical interaction, 
when extended to the Hilbert space of infinite dimension, is proportional to the second virial 
coefficient. 

Very recently, Wu solved the problem of the distribution function for Haldane's frac- 
tional statistics. He found a general form of the equations for the distribution function for 
an arbitrary statistical interaction and discussed the thermodynamics of such a gas. Fur- 
thermore, Bernard and Wu found the explicit form of the statistical interaction in the 
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case of interacting scalar particles in one dimension. In the particular case of bosons they 
showed that, as the amplitude of a local delta-function interaction changes from zero to 
infinity, the distribution function evolves from the Bose-Einstein to the Fermi-Dirac form. 

In this paper we introduce the spin degrees of freedom into the problem and determine 
the statistical properties as well as the statistical interaction for particles in two situations. 
We consider first the Hubbard model in the space of one dimension and with an infinite 
on-site Coulomb repulsion. In this limit, we show rigorously that the charge excitations 
(holons) obey the Fermi-Dirac distribution, whereas the spin excitations (spinous) obey 
the Bose-Einstein distribution. The boson part leads to the correct entropy (/cBln2 per 
carrier) in the Mott insulating limit. As a second example, we express the statistical spin 
liquid partition function with the help of the statistical interaction concept. These two 
examples represent nontrivial generalizations of Haldane's fractional statistics to particles 
with internal symmetry such as spin. In both cases an explicit form of the multicomponent 
statistical interaction is required. We show that the statistical distributions evolve to some 
novel functions when the interaction between the particles diverges. For the spin liquid 
case the form of the distribution functions are also presented for intermediate values of the 
dynamical interaction. In both cases the nonstandard statistics is due to the interaction 
between the particles. 

II. STATISTICAL INTERACTION FOR THE HUBBARD MODEL IN ONE 

DIMENSION 



We consider first the one-dimensional system of particles with a contact interaction. One 
of the simplest models of interacting spin one-half particles was introduced by Hubbard . 
The Hamiltonian in this case is 



A. Thermodynamic limit for the Bethe ansatz equations (U 



oo) 




(1) 
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where t is the hopping integral between the nearest neighboring pairs < i,j > of lattice 
sites, and U is the on-site Coulomb repulsion when the two particles with spin up and down 
meet on the same lattice site. We set t = 1. This model was solved in one dimension by 
Lieb and Wu 0. The solution is given by the set of the Bethe-ansatz equations determining 
the rapidities {ki}, {Aq,}; i.e., 

2tt 1 

—Ij = kj - -Y,'S){2 sin kj - 2A^) 

(2) 

271 1 ^ 1 ^'^ 

— Ja = Aa - 7 E 0(2Aa - 2 sin A;^) + 7 ^ e(A„ - A^g) - ^ A^(5a,,a^, 

^ ^ j=l ^ 13=1 (3=1 

where N is the total number of particles in the system, M is the number of particles with 
spin down, L is the length of the chain, j = 1,...,N, and a = 1,...,M. Ij is an integer 
(half-odd integer) for M even (odd), and Ja is an integer (half-odd integer) for N — M odd 
(even). The phase shift function Q{p) is defined by 

e(p) = -2tan-i(|). (3) 

The second set of equations (2) was written in the form better suited to our purposes. The 
basis in the Hilbert space which diagonalizes the Hamiltonian (1) is called the holon-spinon 
representation. 

The Bethe-ansatz equations can be rewritten in such a way that all dynamical interac- 
tions are transmuted into the statistical interaction We determine explicitly the sta- 
tistical interaction in the case of the Hubbard model. Our method is a straightforward 
generalization of the Bernard and Wu result |^ and is valid in the case of infinite interaction 
only. In this limit, the charge and the spin excitations decouple and there are no bound 
states in the system UTO). In other words, all bound states in the upper Hubbard subband 



are pushed out from the physical many-particle Hilbert space. 



In the large U limit the Bethe ansatz equations read [11 
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M 



27r 1 
-/, = fc, + -E0(2A/.) 

^ ^ /3=1 



M 



M 



(4) 



^ /3=1 /3=1 

We rewrite these equations in the thermodynamic hmit, i.e. for —>■ oo, L —>■ oo, and 
N/L = const. We divide the range of the momentum k's and A's into the intervals with an 
equal size Ak and AA, as well as label each interval by its midpoints ki and Aq,, respectively. 
We treat the particles with the momenta in the i-th or the a-th interval as belonging to the 
i-th or the a-th group. As usual, the number of available bare single particle states are = 
LAk/2TT and = LAA/27r. These numbers follow from the decomposition of the Bethe- 
ansatz wave function in the U —>■ oo limit Next, we introduce the distribution functions 
(the densities of states) for the roots kj and A^ of the Bethe-ansatz equations (4). Namely, 
we define Lp{ki)Ak = as the number of k values in the interval [ki — Ak/2, ki + Ak/2], and 
La{kp)AK = as the number of A values in the interval [A^ - AA/2, A/3 + AA/2]. Hence, 



the two quantities 2T[p{ki) = N^/G'^ 



and 27rcr(A/3) = N^/G^ = rip, are respectively the 



occupation-number distributions for the holon and the spinon excitations in the Hubbard 
chain. In effect, the Bethe ansatz equations in the intervals Ak and A A take the form 



^I{k,) = fc, + E e(2A/3)cr(A/3)AA 
^ J(A,) = A„ - ^e(2A,) + E 0(Aa - A/3)a(A/3)AA 

L L p 

-LEA/35A,,A,a(A/3)AA. 

The function p{k) does not appear explicitly in the large U limit. 

The numbers of accessible states in each of the i-th and the a-th groups are 



(5) 



(6) 



D'^{{Nn, {A^|}) = JiAa + AA/2) - J(A« - AA/2). 
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(7) 



Using the continuous form (5) of the Bethe ansatz equations we find that = Lp^{ki)Ak 
and = Lpl{Aa)AA, where in the thermodynamic hmit {Ak — > , AA — > 0) we have 
respectively the total densities of states for charge and spin excitations 

Ptik) = ^, (8) 



and 12 



"^^^ ' 2n 2nL dA 2-k J ^ ' dA J ^ ' dA ^ ' 

Substituting the form (3) for 0(p) to the derivative dQ/dA one can easily find that in the 
U —>■ oo limit 

p?(A) = i- + a(A). (10) 



To derive (10) we utilized the fact that cr(A) is a flat function of A in the large U limit [|T3 
We see that the numbers of accessible states for the holons and the spinous in the U —>■ oo 
limit are independent of each other. This result, as we show in the following, leads to the 
decomposition of the partition function into the holon and the spinon parts. 



B. Statistical interaction for the Hubbard chain 

We define the statistical interaction and the total number of states for spin one-half 
particles. For that purpose we work in the basis in which the Hamiltonian is diagonal, i.e., 
we choose the holon-spinon representation and observe that the dimensions Df and of 
the one-particle Hilbert spaces for the particle in the i-th or the a-th groups are functional 
of both {N^} and {N^}, i.e. = D^{{N^}, {Nl}), and = D^({iVf }, {ATJ}). Namely, 
starting from the Haldane definition ||^ of the change of the number of the accessible states 
and adopting it to the present situation we obtain 

m = - E 9--^N^ - E aZ^K, (11) 
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= - E daA^j - E 9':,^Nh (12) 

3 13 

where the four g parameters are called the statistical interactions. These difference 
equations can be transformed to the following differential forms: 

(-^?;)-^^+«)-^^ = 2, (13) 

0^=2. (14) 

This set of equations establishes the generalization of Haldane's equations for the number 
of accessible orbitals of the species a in the case of particles without internal symmetry. As 
before |@, statistical interactions {g} do not depend on the occupations A*"^ and N^, since 
otherwise the thermodynamic limit would not be well defined. The factors 2 in the r.h.s. of 
(13) and (14) are irrelevant because they can be incorporated into the g parameters. Then 
the solutions of equations (13) and (14) are 

mim. {N-,}) = G° -j:9--n; -y.9zk. (is) 

j a 

Dim), {iV|}) = Gl -Y.9IW3 -Y.9'.W,- (16) 

i /9 

One should note that these solutions are well defined also in the boson limit, since then the 
corresponding g parameter(s) vanish. The relations ({0}, {0}) = and D^({0}, {0}) = 
G° express the boundary conditions for this problem; the values and Gl represent the 
maximal values of available one-particle states in the situation when the holon and the 
spinon bands are empty. 

Additionally, the total number of microscopic configurations with the numbers {AJ} and 
{Np} of holon and spinon excitations is given by 

mm - 1)! iii mm - 1)! ' ^ ^ 

In this expression the two products are in general interconnected via the relations (15) and 
(16). Each of the factors is defined in the same manner as in Ref. 0. In the fermionic 



bookkeeping for Ij and Ja the same fl is obtained with the number of accessible states in 
the i-th and a-th groups taken to be 



j a 

(19) 

Rewriting these equations for each of the intervals Ak and AA we easily find that in the 
Ak and AA limits these four types of statistical interactions reduce to the following 
form 

g''%k,k')=5{k~k'), (20) 



g'=%k, A) = g^^iA, k) = /^(A, A') = 0. (21) 

Thus, the vanishing g functions in (21) simplify the expression (17) for the total number of 
available configurations, which is then 

The statistical weight Vl factorizes into the holon (fi^) and the spinon (fi") parts. This, 
once again, expresses the fact that the spin and the charge degrees of freedom decouple in 



the U ^ oo limit ||T0|- As a consequence, the entropy of the system is a sum of the two 
parts S = S"^ + = kslTiQ'^ + ksiTiQ^ , where the corresponding expressions calculated per 
particle are 



1 ^ 



= -kB^ E Klnri- + (1 - <) ln(l - <)] , (23) 



and 



1=1 



M 



= -kBjT- E K Inn^ - (1 + <) ln(l + <)] , (24) 

^^a a=l 



where Na is the number of atomic sites. 

We recognize immediately that the holon contribution to the system entropy coincides 
with that for spinless fermions, whereas the spinon contribution reduces to locahzed-spin 
moments (A;Bln2) in the Mott-insulator hmit and in the spin disordered phase, i.e. when 

= 1 and M = Na/2. In general, one may say that S*^ provides the entropy of charge 
excitations (and vanishes in the Mott insulating limit n'^ = 1), whereas represent the spin 
part of the excitation spectrum. This demonstrates again that the holon (charge) excitations 
are fermions and the spinon (spin) excitations are bosons. In the U oo limit considered 
here the Heisenberg coupling constant (J = At'^/U) vanishes and the spin wave excitations 
do not interact with each other In other words, they are dispersionless bosons. Also, 



the charge excitations are spinless fermions. The total entropy of the system reduces in the 
Mott insulating spin-disordered limit to S = A;5ln2. This value is different from that for 
the Fermi liquid in the high-temperature limit, which is 2/csln2. This difference confirms 
on a statistical ground the inapplicability of the Fermi liquid concept to the Hubbard model 
in one dimension. 

III. STATISTICAL INTERACTION FOR THE SPIN LIQUID 

In this section we derive the statistical interaction for the so-called statistical spin liquid. 
This concept was introduced in to describe the thermodynamic properties of strongly 
interacting electrons. The basic assumption in this approach is to exclude the doubly oc- 
cupied configurations of electrons with spin up and down not only in real space but also in 
reciprocal space (with given k). This assumption leads to a novel class of universality for 
electron liquids. Its thermodynamics in the normal, magnetic and superconducting states 
were examined in the series of papers 0,0. A justification of this approach has as its origin 
in the concept of the singularity in the forward scattering amplitude due to interparticle 
interactions. Namely, it was noted by Anderson |]T6| and by Kveshchenko [|17| that in two 



spatial dimensions this amplitude may be divergent either due to the Hubbard on-site repul- 
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sion, or due to the long-distance current-current interaction mediated by the transfer gauge 
fields. With the assumption that in those situations the wave vector is still a good quantum 
number, one can write down the phenomenological Hamiltonian describing such liquid in 
the form 

H = E(ek - (rh)N^^ + f/, ^ N^^N^^. (25) 

kcr k 

In this model is the dispersion relation for the particles with the wave vector k moving 
in the applied external magnetic field h, A^ko- is the number of electrons in the state |kcr >, 
and cr = ±1 is the projected spin direction. The number of double occupancies in given k 
state is A^kd = ^kt^kj,- The nonvanishing A^kd causes an increase of the system energy by 
Us > for each doubly occupied k state. Finally, we will put Ug oo because this model 
is to represent the situation with the singular forward scattering amplitude. It turns out 



that this exactly solvable model belongs to the class of models with Haldane's fractional 
statistics, as shown below. In this case the statistical interaction is proportional to 5kk' in 
k space, but is a nondiagonal matrix in the extended spin space. 

To prove this we define the total size of the Hilbert space of the many-particle states 
determined by the number of physically inequivalent configurations 

V (iVkt - iVkd)!(DkT - 1)! (iVki - N^MDkl - 1)! (iVkd)!pkd - 1)! ' ^ ^ 

Due to the local nature of the interaction in k space we must treat separately the singly 
occupied states as distinct from those with double occupancy in reciprocal space. Then, 
the statistical weight Q expresses the possible ways of distributing A^k^ — A^kd quasiparticles 
over Dko- states and A^kd quasiparticles over the Dkd states. In general, the dimension of the 
one particle Hilbert space for the singly (-Dko-) and the double occupied (-Dkd) states is the 
function of the number of other quasiparticles {A^ka} and {N^a} 0Jl^ i.e., 

DU{N^.p}) = Gl -Y,g^^p{k,k'){N^>^ - 5^^^ (27) 

where a and /3 label the configurations j, these states define the extended spin space. 
Note that in difference with Eqs.(15) and (16) we define here the boundary conditions via 
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the relations -Dka({A^k/3 = <5a/3<5kk'}) = G^, i.e. the maximal dimension of the single-particle 
Hilbert space is defined for an occupied configuration in each category, not for an empty 
one. This new conditions are equivalent to the form appearing in Eqs. (15) and (16), in the 
thermodynamic limit. Since, the Hamiltonian (25) does not mix different momenta, we find 
the general solution for 5'a/3(k, k') in the form 



^a/3(k,k') = (5kk' 'S)g, 



(28) 



Hence, the statistical interaction is diagonal in k space for this model. Next, to get the 
exact solution of the Hamiltonian (25) we choose 

/ \ 

1 1 -1 



1 



(29) 



and substitute (29) and (28) into (27) we find that the total size of the many-particle Hilbert 
space is given by 



(30) 



This result is exactly the same as that obtained in Ref. [0 (cf. Appendix B). Therefore, we 
conclude that this model also belongs to the class of models with Haldane's statistics. In this 
case, the change in the distribution functions are not due to the phase shift between different 
momenta but rather due to the mutual (dynamic) interaction between quasiparticles with 
the same k but different spin. The interaction pushes some of the states upward in energy, 
leading to the following form of momentum distribution functions 



Gl 



e^f/.e/3(^k-M) cosh(/3/i) 
1 + e^f^=e/3(^k-M)[e/3(^k-M) + 2 cosh(/3/t)] 



[1 + atanh(/5/i)]. 



Gi 1 + e^^^ e/^(^k-/^) [e/5(^k-M) + 2 cosh(/3/i)] ' 



(31) 



(32) 
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which are easily obtained by minimizing the thermodynamic potential with respect to A^k^ 
and A^kd separately [|1^. It is easy to show that those distributions evolve from the Fermi- 
Dirac function to the statistical spin liquid distribution when the Us changes from zero to 
infinity §. 

The limit Ug ^ oo represents the physical situation in which A^kd = 0. In other words, 
there are no double occupancies in k space. All states are singly occupied by the quasipar- 
ticles with either spin up or down, or empty. In this limit, the statistical interaction (28) 
reduces to the 2x2 matrix form 



1 1 



(33) 



Then, the statistical weight is 



(34) 



Such liquid is called the statistical spin liquid. This is a novel class of quantum liquids, 
similar, in some respect, to the Bethe-Luttinger liquid discussed above. For example, because 
of the mutual interaction between spin up and down particles, half of the total number of 
states (2A'a) are pushed out of the physical space in the Us ^ oo limit. Therefore, the 
entropy of the statistical spin liquid in the high-temperature limit is the same as in the case 
of the Hubbard chain with the infinite interaction because the entropy of the system in this 
temperature limit is given in terms of the degeneracy of the state only [|2l[. In particular, 
the entropy in the statistical spin liquid for N = Na equals ksN ln2, which is the correct 
value for a Mott insulator. Also, the high-temperature value of a thermopower is the same 
for both liquids It was also shown |19| that the magnetization of statistical spin liquid 
has the same form as that of the Hubbard chain with the infinite repulsion, i.e. that for 



localized moments |]TT| . However, the direct consequence of the statistical interaction is also 
a breakdown of the Luttinger theorem: the volume enclosed by the Fermi surface is twice 
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that for the Fermi hquid. This arises because of the differences in the microscopic character 
of the single-particle excitations in these two liquids. 

IV. CONCLUSIONS 

In this paper we considered statistical properties of the two model system: The Hubbard 
chain with infinite repulsion and the statistical spin liquid. We determined the form of 
the Haldane statistical interaction in each case. In the one-dimensional Hubbard model a 
novel distribution function emerges due to the presence of the phase shift between pairs 
of states with different rapidities {Aq,} and {fcj}. In the U ^ oo limit, when all bound 
states are excluded, the charge excitations (holons) behave statistically as fermions, and the 
spin excitations (spinous) behave as bosons. The holon have a simple energy dispersion 
coinciding with the bare band energy, whereas the spinous are dispersionless. An equivalent 
alternative approach |2y] is based on the fermionic representation of both the holon and the 
spinon degrees of freedom; in that approach the spinous and the holons acquire complicated 
forms of the effective dispersion relation. The latter approach allows for a generalization 
of the treatment of the Hubbard chain for finite U, the limit unavailable with our present 
analysis. Nonetheless, our analysis is also applicable to other models when the separation 
into the charge and the spin degrees of freedom occurs (cf. the Kondo problem |p2|| ). 

In the statistical spin liquid case the mutual interaction between quasiparticles with the 
same k leads to the exclusion of the double occupied configurations in reciprocal space. In 
that case the statistical interaction is diagonal in k but has a nondiagonal structure in spin 
space. In that case the statistical distribution changes with growing interaction from the 
Fermi-Dirac form to the spin liquid form [0. 

It is interesting that those two models of particles with the internal symmetry can be 
classified as the models with the fractional statistics in the Haldane sense. By contrast to 
the case with scalar particles 0, the distribution functions in the present situation take 
the form of either holon-spinon or the spin liquid distributions. Those possibilities arise 
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only when the particles have some internal symmetry (spin,colour). One may also say that 
in those cases the statistical interactions have a tensorial character in space in which the 
Hamiltonian is diagonal. 
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